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Введение. Статья посвящена вопросу анализа чувстви- 
тельности функционала качества воды к параметрам мо- 
дели искусственной водной экосистемы биологической 
очистки сточных вод, сформулированной в виде начально- 
краевой задачи для системы уравнений реакции- 
конвекции-диффузии, описывающей динамику биогенов. 
Использование многомерных нестационарных моделей 
водных экосистем с подробным описанием биохимиче- 
ских реакций осложняется большим числом параметров, 
значение которых необходимо получать эксперименталь- 
но. Для упрощения процесса обеспечения данных была 
исследована чувствительность функционала качества 
очистки к параметрам модели. 

Материалы и методы. В работе проведён краткий обзор 
существующих методов оценки, на основе методов теории 
малых возмущений и сопряжённых уравнений сформули- 
рован алгоритм исследования чувствительности функцио- 
нала, характеризующего качество воды, к параметрам 
модели водной экосистемы. 

Результаты исследования. Получен алгоритм исследова- 
ния чувствительности функционала к параметрам модели 
водной экосистемы. Приведён практический пример при- 
менения рассмотренного метода для задачи оптимального 
управления аэротенком-вытеснителем — элементом 
очистного сооружения для биологической очистки сточ- 
ных ВОД. 

Обсуждение и заключения. На основе проведённого ана- 
лиза выявлены наиболее и наименее значимые параметры 
математической модели искусственной водной экосисте- 
мы аэротенка-вытеснителя, являющейся составной частью 
задачи оптимального управления аэрацией. Применение 
разработанного алгоритма оценки чувствительности воз- 


можно и для других процессов, связанных с массоперено- 
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сом реагирующих веществ, таких как решение обратных 
источниковых задач динамики и кинетики газовых приме- 
сей и аэрозолей в атмосфере, моделирование биологиче- 
ских процессов в живых организмах, управление процес- 
сами массопереноса в аппаратах химических технологий и 
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Введение. Рассматривается задача оптимального управления процессом биологической очистки сточ- 
ных вод на основе математической модели, описывающей динамику концентраций биогенов в очистном со- 
оружении. При решении задач управления экологическими системами на основе сложных многомерных и не- 
стационарных моделей с подробным описанием биохимических реакций возникают сложности, связанные с 
большим числом параметров, значение которых необходимо получать экспериментально [1]. Для упрощения 
процесса обеспечения входных данных многопараметрической модели авторами исследуется чувствительность вы- 
бранной значимой скалярной характеристики (функционала качества очистки сточных вод, формирующего ограниче- 
ние на переменную состояния в задаче оптимального управления) к параметрам модели. Путём этого исследования 
параметры можно разделить на наиболее значимые, для определения которых рекомендовано проведение натурных 
экспериментов, и наименее значимые, значения которых могут быть взяты приближённо. 

Наиболее простым способом оценки чувствительности функционала к изменениям параметров много- 
параметрической модели является метод стохастического анализа (как линейного, так и нелинейного) [2-4]. 
Кроме этого, возможно применить метод, основанный на теории малых возмущений и сопряжённых уравнений 
[2, 5—7]. С вычислительной точки зрения этот алгоритм имеет преимущества по сравнению с прямыми метода- 
ми, так как вместо многократного решения уравнений для множества параметров модели при расчёте компо- 
нентов градиента функционала предполагается одноразовое решение прямой и сопряжённой задач и вычисле- 
ние интегралов, характеризующих оценку влияния изменения параметров к функционалу. 

Данный подход реализован авторами в настоящей работе. Применение алгоритма оценки чувствительно- 
сти, сформулированного ниже, возможно для любых процессов, связанных с массопереносом реагирующих ве- 
ществ, таких как решение обратных источниковых задач динамики и кинетики газовых примесей и аэрозолей в 
атмосфере, моделирование биологических процессов в живых организмах, управление процессами массопереноса 
в аппаратах химических технологий и т. д. Возможно применять данный подход и для других классов задач. 

Описание метода и алгоритма. Рассматривается задача, описывающая распространение и массоперенос 
неконсервативных примесей, представляющая собой краевую задачу для полулинейной системы уравнений реак- 
ции-конвекции-диффузии: 
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(х,х,,х,) <> (х,у,2) — пространственные координаты, К+А, — сумма коэффициентов молекулярной и турбу- 


лентной диффузии (в направлении х, ) соответственно. 
Пусть имеется некоторый интересующий нас функционал: 
Л = ( р,э ) | 
где знак <, ) обозначает скалярное произведение в области Тх© (Т: {1=[щ, Т}, Ф — пространственная расчётная 


область), т. е. (",й) = | "(х)й(х) аа! ‚ здесь 7(х), й(х) — некоторые произвольные функции, заданные в обла- 
Тхо 


сти 1х0; р={р,} — вектор функция, 1-й элемент которой равен 1 в значимых по веществу 7 зонах 62, и равен 0 
в остальной области, т. е.: 
| при хеО, ХТ, 
в, (х) = (2) 
О при хЕО, ХТ. 


Вдали от точек бифуркации решения системы (1) малые изменения параметров порождают возмущён- 
ную задачу, решение которой можно представить в виде ряда: 


_ 25[1.^ 
5(.^) = 5 »+ 5 _^5А 
2 2 ие ал р) 
где 0) — вариации параметров, 9%. =^,—^, = {0 @%. 
051 


Функция чувствительности по [-му параметру ‚ является решением системы уравнений модели 
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(1), линеаризованной относительно невозмущённой траектории 5 =5 (2. | 
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задаёт систему уравнений линейного приближения, где запись |. обозначает матрицы частных производных 
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правой части системы. Первая вариация функционала определяется по формуле: 
&Л, =(55,р)=(в,60), 1=1,Ё, (4) 
где « — решение задачи с оператором А”, сопряжённым к 4. 

Определение вариации функционала (4) в рамках прямого подхода предполагает расчёт функций чув- 
ствительности, а для многопараметрической модели — многократное решение начальной задачи для системы 
уравнений (3). Использование решения сопряжённой задачи А” позволяет вместо многократного решения урав- 
нений линейного приближения для множества параметров основной задачи при расчёте компонентов градиента 
функционала .Л реализовать подход, который предполагает одноразовое решение прямой и сопряжённой задачи 


и нахождения скалярного произведения: 


7 _ ©, се ‚5 1=12. (5) 
07. д. 





Алгоритм оценки чувствительности функционала к изменениям параметров моделей при помощи тео- 
рии малых возмущений и сопряжённых уравнений предполагает выполнение следующих этапов: 
1. линеаризация полулинейной системы уравнений реакции-конвекции-диффузии (1); 
формулировка задачи, сопряжённой к линеаризованной; 


нахождение решения прямых и сопряжённых задач, 


| 
вычисление скалярного произведения (5) для каждого параметра /, ‚, [=1Г; 


А 


оценка полученных компонент градиента функционала. Чем больше соответствующая компонента, 


[ 
тем большее влияние оказывает изменение параметра Л на функционал. 
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Постановка задачи управления биологической очисткой сточных вод и результаты численных 
экспериментов. Определим влияние вариации различных параметров модели экосистемы аэротенка- 
вытеснителя на функционал, характеризующий эффективность процесса очистки сточных вод от загрязняющих 
веществ. Под аэротенком-вытеснителем здесь понимается элемент очистного сооружения, в котором воды, за- 
грязнённые аммонийным азотом (№), являющимся наиболее опасным для организма человека веществом, сме- 
шиваются с активным илом (А), а также аэрируются кислородом (ОХ), необходимым для протекания биохими- 
ческих реакций. В ходе прохождения «смеси» по аэротенку происходит очищение сточных вод от загрязняю- 
щего вещества в процессе наращивания биомассы активного ила [8, 9]. Схема аэротенка-вытеснителя изобра- 
жена на рис. 1. 





кислород 


Рис. 1. Схема аэротенка-вытеснителя 
Ето. [. эсйете оГаегайоп апк 


Сформулируем задачу оптимального управления аэрацией элемента очистного сооружения. 
Управлением является мощность источника кислорода Роу (+) . Функция цели задачи обеспечивает ми- 


нимизацию затрат на поступление кислорода в водную систему аэротенка-вытеснителя в течении периода вре- 
мени Т: 


ЕЕ, (Е) 4 —шш, 


где &, — затраты, необходимые для подачи кислорода объёмом Ру (+) на единицу площади источника в мо- 





мент времени [=&,Г. 
Переменными состояния задачи являются концентрации №, А и ОЛ: 5, 5., эох. Оптимизационная 
модель включает в себя следующее ограничение на переменную состояния: 
(р» (х),5» -5»рк }<0, хеО, (6) 
где бу ррк — предельно допустимая концентрация аммонийного азота; р» (х) — функция типа (2), характери- 


зующая значимую область (соответствует зонам водозабора из сооружения). 
Управление удовлетворяет ограничению на знак: 
ВОГ 
Решение задачи оптимизации должно удовлетворять краевой задаче, описывающей массоперенос и 
трансформацию неконсервативных примесей (уравнениям состояния): 











я +4 (и-5.)-2$, И с 2 и А 
К +54 ки о (7) 
р о. Е.Г, 280, 
05 - 5 т 
и +4 (и-5„}-Абу =655, — и 2 баз Е=А,Г, хе, (8) 
Ку+5+—- 
№ 
05 - 5 т -- 
ое (и. ох бок = 5, ГЕК, Г, хЕ@); (9) 
О Кох +5ох 
0, 1={А, №, ОХ}, хе ОО, ее Хх ох; (10) 
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ток = 5 або) ФЕТ: МЕХ. (11) 

9=57, 1={А,М, ОХ}, хЕО, 1=4. (12) 

В вышеприведенных уравнениях |, — максимальная скорость роста активного ила; К,, {={А, МОХ} — кон- 
станты полунасышения; К", {={А, №, ОХ} — константы ингибирования роста активного ила; ф, — коэффи- 
циент внутривидовой конкуренции; /, — скорость отмирания; 5) — скорость выделения аммонийного азота в 
процессе жизнедеятельности активного ила; К, — коэффициент, равный отношению скорости растворения 


кислорода к коэффициенту диффузии в направлении нормали №. Задача состоит в оценке чувствительности 
функционала из ограничения (6) к параметрам модели. 
Для реализации изложенного выше алгоритма необходимо сформулировать задачу нахождения вектора 


(1) = (5. (#),5» (1) 5х (1) ‚ то есть сопряжённую к краевой задаче [10], являющейся линейным приближени- 


ем к задаче (7)—(12). 
Сопряжённая задача имеет вид: 


* 
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М М№ т 
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——> 4 бы +45, — 
———[. 25 
ИИ а 
ИиЗох А Ку . 
— т р-он (14) 
а в” —_ — 
—__ —_ 
КУ ки (Кох +8ох } Кб (Кох + Зих } 
АУ К» ы 
о м = Рь, 
о В Ким 
ыы М 
о а р 
-—1: 4% бох +Абох — 
5,5, а. р 
__ Ре А __ Рим ОХ -А 5. - (15) 


К ++ (Кох + ох } К (Кок +5 


М№ 


5 о 
_| ___ Им Ау Ит ох 5 5, =0; 


(Кох а бох (Кох + а. 
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их, + > =0, 1={ А, №,ОХ1, хебО, 1=&,Т, хеХ ох; (16) 
9 801. МОХ, ЕО ЕЕЕ. (17) 
В системе (13)-—(17) 5, ‚5 1= { А, МОХ } — среднее значение концентрации, полученное при решении 


системы ('7)—(12) с невозмущёнными параметрами. 
Далее определим чувствительность функционала из ограничения (6) к изменениям следующих пара- 


метров /,: К,, [= {№ ‚ОХ } о - [ д; Фа; М» - Результаты вычисления скалярных произведений (5) сведены в 


таблицу 1. 
Таблица 1 
ТаБЕ 1 


Чувствительность функционала эффективности процесса биологической очистки 


Епслепсу ипсНпопа! зепяй»йу ор оозса![ ритйсаНоп 





Обсуждение и заключения. Анализ чувствительности функционала показал, что параметры, характе- 
ризующие константу полунасыщения по кислороду (Ко, ), константу полунасыщения по аммонийному азоту ( 


К»), максимальную скорость роста активного ила (и„), скорость отмирания активного ила (/,) и коэффициент 
внутривидовой конкуренции (ф,) оказывают наибольшее влияние на вариацию функционала эффективности 


очистки сточных вод и их значения следует определять экспериментально. В тоже время изменения парамет- 
ров, характеризующих константу ингибирования роста активного ила ( К’) и скорость выделения аммонийного азо- 


та в процессе жизнедеятельности активного ила ( 6%, ) практически не влияют на изменение функционала, характе- 


ризующего эффективность очистки сточных вод, и, следовательно, при расчетах можно использовать средние 
значения этих параметров. 
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